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We calculate Wilson loops of various sizes up to loop order n = 20 for lattice sizes of L 4 (L = 
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plaquette data and our perturbative result we estimate a value of the gluon condensate (^GG). 
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1. Introduction 

Since the introduction of the non-perturbative gluon condensate by Shifman, Vainshtein and 
Zakharov [jl|] there have been many attempts to obtain reliable numerical results for this quantity. 
Soon it became clear that lattice gauge theory provides a promising tool to calculate it from Wilson 
loops. In [Q] the plaquette was used whereas larger Wilson loops have been investigated in [^]. 
From the plaquette P the non-perturbative gluon condensate (f GG) is conventionally derived from 
the relation 



MC — r pert 



Pnprt Q> 



K 



36 



-bo; 



Pig) 



(-GG), (1.1) 
n 



where bo is the first coefficient of the /3 -function and P MC is the plaquette measured in Monte Carlo. 
In it is assumed that the non-perturbative part scales like the fourth power of the lattice spacing 
a. However, there were speculations that there could be non-perturbative contributions which scale 
like a 2 |@]. In the last decade the application of Numerical Stochastic Perturbation Theory (NSPT) 
[^] pushed the perturbative order of P pert up to order n = 10 [§] and even n = 16 [0]. This strongly 
supports to use ( |1.1[ ) for the determination of (^GG). 

Besides the determination of {^GG} there is a general interest in the behavior of perturbative 
series in QCD (for a recent investigation see [||]). Observable quantities can be written as series of 
the form 

<2~£a„A'\ (1-2) 

where X denotes some coupling. It is generally believed that these series are asymptotic, and 
assumed that for large n the leading growth of the coefficients a„ can be parametrized as ^ 

fl»~C 1 (C 2 )"r(n + C 3 ), (1.3) 

i.e., they show a factorial behavior. Using the technique of NSPT one reaches orders of the per- 
turbative series where a possible set-in of this assumed behavior can be tested. There is a recent 



paper of Narison and Zakharov [10] where the authors discuss the difference between short and 
long perturbative series and its impact on the determination of (^GG). 

In this paper we present perturbative calculations in NSPT up to order n = 20 for Wilson loops 
for lattice sizes L 4 with L = 4, . . . , 12. The computation for L = 12 were performed on a NEC SX-9 
computer of RCNP at Osaka University, all others on Linux/HP - clusters at Leipzig University. 
We calculate the Wilson loops in quenched QCD with plaquette gauge action. 



2. NSPT calculation up to n = 20 

NSPT allows perturbative calculations on a lattice up to loop order n which never will be 
reached by the standard diagrammatic approach. The algorithm is introduced and discussed in 
detail in [II]] - we will not present it in this paper. We only want to point to some essential 
topics: 

• The computer implementation of NSPT requires the discretization of the so-called (rescaled) 
Langevin time % 

z^T k = ke/ji, k = 0,l,2,.... 
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(g = 6/f3 is the bare lattice coupling). Practically, this means that the corresponding quan- 
tities are measured for different small but finite e. The final result is obtained in the limit 
£ — > 0. This must be done with great care in order to obtain reliable numeric results. 



The connection to infinite volume is achieved by the limit L 
extrapolation of the corresponding finite L results. 



• oo which requires an additional 
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Figure 1: Extrapolation e — > for W\\ for 1-loop (left) and 20-loop (right) forL = 8. 

In Fig. [I] we show the extrapolation e — > for lattice size L = 8 for a plaquette, where we use 
a general quadratic ansatz in e for the fitting function. 

We write the general expansion of a Wilson loop of size N x M in terms of the bare lattice 
coupling g as 

20 / ^ 

W nm =Y, w nmS ■ (2-1) 



Depending on the loop-size (N,M) we found alternating signs for the perturbative coefficients wj^ 
for smaller n whereas for larger n they turn into a smooth asymptotic behavior. An example is given 
in Fig. | (left) for L= 12. 



A typical extrapolation to L — > oo for the plaquette is shown on the right side of Fig. g. Bali [ |12| ] 
has computed one- and two-loop contributions to Wilson loops of various sizes in the standard 
diagrammatic approach for finite L. A comparison of our one- and two-loop NSPT results with his 
results is given in Table [IJ Based on the results given by Bali we fixed the functional dependence of 
the L — > oo extrapolation. However, it should be empasized that this extrapolation becomes worse 
for larger loop sizes (N,M). 



3. Perturbative series at large order 

The order of perturbation theory we have reached in our calculations allows to study the large 
order behavior and to test some models concerning the n— dependence of the coefficients. This is 
essential in order to compute the perturbative part of the Wilson loops as precise as possible. In 
order not to interfere with possible extrapolation (L — > oo) effects we investigate this for finite L. 
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Figure 2: Coefficients for various Wilson loops (left). Extrapolation L — > °o for « = 10 (right). 
Table 1: Comparison of one- and two-loop results for NSPT and standard approach 
W NN L NSPT (1-loop) Bali(l-loop) NSPT (2-loop) Bali (2-loop) 



W22 


4 


-0.87468(13) 


-0.87500 


0.10404(07) 


0.10406 




6 


-0.90752(12) 


-0.90762 


0.11830(10) 


0.11837 




8 


-0.91164(08) 


-0.91141 


0.12008(08) 


0.11993 




12 


-0.91259(03) 


-0.91261 


0.12038(04) 


0.12038 


W33 


6 


-1.50088(30) 


-1.50093 


0.60906(34) 


0.60866 




8 


-1.52873(23) 


-1.52803 


0.63693(23) 


0.63632 




12 


-1.53526(47) 


-1.53533 


0.64370(13) 


0.64360 




8 


-2.14128(44) 


-2.14016 


1.52351(70) 


1.52331 




12 


-2.16950(24) 


-2.16922 


1.57178(60) 


1.57006 



3.1 Heuristic model 



In [13] the authors propose to use a series expansion for a quantity which shows a power-like 
singularity 



r(n - 7) 



w u ~ a-us 2 y=y — - 



lug 



,2\n 



2n 



(3.1) 



From ( J3- 1[ ) one derives the ratio of successive coefficients c n as (slightly modified by a parameter s 
to account for a small curvature) 



r n = c n /c n -\ = u\\ 



1 + 7 



(3.2) 



4 



Wilson loops at very high order LPT 



H. Perlt 



In a Domb-Sykes plot - r n plotted against l/n - this is almost a straight line. In Fig. [3] one observes 
that r n for W\ \ follows this simple functional form almost ideally. 

However, the corresponding curves for larger Wilson loops of moderate size have a more 
pronounced non-linear dependence on l/n as can be seen in Fig. ||| This suggests to generalize 
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Figure 3: Domb-Sykes plot for various Wnm together with their fits (3.3). 



ansatz ( p.2| ) by adding an extra power in n (for a detailed discussion see Q14[]) 



n 2 + (s-q-\)n + t 

c n /c n -]_=u — — . (3.3) 

n(n + s) 



For t = relation ( p.3| ) is identical to (p.2|). It gives a hyperbola in a Domb-Sykes plot. In this paper 
we assume that the intercept u has a universal value for all loop sizes (N,M). It is determined from 
Wn which has been computed most precisely. The other parameters (q,s,t) depend on (N,M). The 
corresponding curves are shown in Fig. |5[ They are obtained from the fit ansatz ( [3.3[ ) where the 
parameters are determined in the interval 5 < n < 20. In this region the perturbative coefficients 
of the considered Wilson loops show a common asymptotic behaviour as can be seen in Figure 
| (left). 

There were speculations that already at order n = 10 the perturbative coefficients show a fac- 
torial growth due to renormalon contributions [01 0] (for a detailed investigation of this point see 



also [||]). For the plaquette we plot in Fig. |] the ratio r„ over n for the ansatz (3^) (HRS) and the 



renormalon inspired model as given in [Q, |6|] (BDMO). We do not observe a factorial growth, at 
least in the region n < 20 and for our lattice sizes. 

3.2 Boosted perturbation theory 

It is well-known that the bare lattice coupling g is a bad expansion parameter due to lattice 
artefacts like tadpoles. There is a hope that by redefining the coupling g into a boosted coupling gh 
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Figure 4: Comparison of r n of the plaquette W\ \ for HRS and BDMO models 



and the corresponding rearrangement of the series a better convergence behaviour can be achieved. 
For the plaquette P = W\\ we use the replacements 



8l 



pert.b 



Ppert(g,n* 



1 + I WtfV" - P P ertA8l,y) =1 + 1 <iU" , (3-4) 



(") „2n 



n=\ 



n=\ 



where n* is the maximal loop order. 

Boosted perturbation theory has been applied to improve the perturbative series for the pla- 
quette for the first time by Rakow [0]. He showed that P pe rt,b(gb,n*) reaches a stable plateau much 
earlier than P per t(g,n*) as a function of n*. Fig. || (left) shows that the boosted coefficients wjf 1 ^ 
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Figure 5: Coefficients for naive and boosted LPT (left). P at j3 = 6.2 as function of n* (right). 

oscillate but rapidly become very small. Of course, one should act with caution in the region of n 
where \W^\ ~ 1(T 7 . The superior convergence behavior for the plaquette is demonstrated in Fig. 
|5] (right) confirming the result in [f7j|. The Monte Carlo result is taken from Jl5| , |TJ]. 
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4. Non-perturbative gluon condensate 



As discussed in the introduction there are speculations whether the difference AP = P, 



pert 



behaves as ~ a or ~ a . We can check this by plotting AP versus a/ro where tq denotes the 
Sommer scale. The functional relation between /3 and ro/a has been taken from Jl7|]. In Fig. || 
AP(a/rr,) is plotted in the infinite volume limit (L — > oo) for both models discussed in the previous 
sections. The MC data points have been taken from [15, [if]]. (The cut-off in the HRS-model data 
for larger a is due to the convergence radius for the coupling determined by the parameter u in 
(3.1).) We make the ansatz AP(a/ro) = C(a/ro) 4 and approximate ( ~fi[g) ) ~ 1- This gives for the 
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Figure 6: AP(a/ro) with fit curves ~ (a/ro) 4 . 

range 0.1 < a/ro < 0.25 

r 4 (-GG) ms = 1.63(9), r 4 {-GG) boosted = 1.80(5). (4.1) 

Fig. H shows that the data are well described by the ansatz ~ (a/ro) 4 over a large range of a. 
Inserting , e.g. ro = 0.5 fm we obtain 

(-GG) H rs = 0.039(2) GeV 4 , {-GG) boosted = 0.043(2) GeV 4 . (4.2) 
K K 

One can try to fit the more general ansatz AP = C (a/ro) 5 to the data. For the boosted model and 
0.1 < a/ro < 0.25 we get 8 = 3.5 ±0.1 which is not too far from 5=4. 

All given errors are purely statistical, some of the systematic uncertainties are at least as large, 



and we are planning a more careful error analysis in the full paper Q14[]. It should be emphasized 
that the determination of the gluon condensate depends less on the assumption of large loop order 
behavior than in earlier investigations where all contributions beyond n = 10 were obtained by 
extrapolation. 

5. Summary 

In this paper we presented the perturbative calculation of Wilson loops of different sizes up 
to loop order n = 20 using NSPT We compared three models to describe the data: a renormalon 
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inspired model (BDMO), a heuristic fit (HRS) and boosted perturbation theory. We found that up 
to order n = 20 the resulting curves show a ~ a 4 behaviour. This supports the claim of Narison and 
Zakharov [10] that a behaviour ~ a 2 is due to perturbative series cut a lower order. The values (| 



for (^GG) found for HRS and boosted PT are larger than obtained in other computations [7| [T3| ]. 
The gluon condensate can also be obtained from larger and/or asymmetric Wilson loops serving as 
an additional check. We hope to come back to this problem in [14]. 
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